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Abstract. We present the study of one year of INTEGRAL data on the neutron star low mass X-ray binary GX 5-1. Thanks 
to the excellent angular resolution and sensitivity of INTEGRAL, we are able to obtain a high quality spectrum of GX 5-1 
from ^5 keV to '--^lOOkeV, for the first time without contamination from the nearby black hole candidate GRS 1758-258 above 
20 keV. During our observations, GX 5-1 is mostly found in the horizontal and normal branch of its hardness intensity diagram. 
A clear hard X-ray emission is observed above ~30 keV which exceeds the exponential cut-off spectrum expected from lower 
energies. This spectral flattening may have the same origin of the hard components observed in other Z sources as it shares the 
property of being characteristic to the horizontal branch. The hard excess is explained by introducing Compton up-scattering of 
soft photons from the neutron star surface due to a thin hot plasma expected in the boundary layer The spectral changes of GX 
5-1 downward along the "Z" pattern in the hardness intensity diagram can be well described in terms of monotonical decrease 
of the neutron star surface temperature. This may be a consequence of the gradual expansion of the boundary layer as the mass 
accretion rate increases. 
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1 . Introduction sources have weaker magnetic fields (about 10^ to 10^ G versus 

10^-10^ G of Z sources), are generally fainter (0.01-0.3 Lfidd 
Low-mass X-ray binaries (LMXRBs) are systems where ^^^^^^ ^^^^^^^ ^^^^-^ ^^^^^^.^ ^.^^^ ^j^^ 2 

^ ; the compact object, either a neutron star (NS) or a black ^j^^p^ ^-^^ ^^^j^^ ^j^^^^ ^^p^^^j z.^o^rces. Besides, 

hole candidate (BHC), accretes matter from a companion ^j^^j^ ^.^.^^ variability properties in a given spectral state are 

with a mass M< IMq. LMXRBs hosting a weakly magne- (jjfferent 

tised neutron star can be bro adly classified in two classes ^, , ^ , ■ , , , . ■ 

rz — : ; — |,„„^ ,^ , The study or these sources in the hard X-ray domain 

(lHasinger & van der Klisll989l) : high luminositv/Z sources and , . • „ • ■ • ■ 

' ; — T : — 71 — ~ f , , , , has proven important especially in the recent years, in the 

low luminosity/AtoU sources. In the colour-colour (CC) and X- ,. , ^ , ■• ^ 7 ■ •, ■• 

. . ,. ,TTTT^^ . „ , liRht or the discovery or hard tails extending up to ^lOOkeV 

ray hardness intensity diagrams (HID), Atoll sources display . , ,,,,t^„ r^.. . . . . 

^ „ ^ , •, in NS LMXRBs. This kind of hard emission was thought 

an upwardly curved branch while Z sources describe an ap- , . ^ t,tt^ . . . , 

f „t,,, , ... , ,• — ~n to be a prerogative or BHCs and had been proposed as 

proximate Z shape. Although two recent studies (MunoetalJ ... . r , r x,tt ; 



miGierlinski fc n;;;;3m have suggested that the clear ifr^P^^^b''"!, ^J^'^TT\''i '' J'' k 

7777^ I .. . : , ' ^ .^^ , iMcChntock & Remillardi2004 . Hard tails discovered in about 

Z/Atoll distinction on the CC diagram is an artifact due to in- 



.„ , , 20 NS LMXRBs of the Atoll class as well as in some Z sources 

complete sampling (Atoll sources, if observed long enough, , , , , . . , . ■ 

.... _ , . . „ showed th at NS binaries as w ell can po wer such an emission 

do exhibit a Z shape as well) many differences remain. Atoll jr — ■ ~, , „. 7^ „ ^„„^, 

jBarret & Vedrenne 1994: Di Salvo & Stella 200 2 . 



Send offprint requests to: Ada.Paizis@obs.unige.ch 5-1 (4U 1758-25) is thought to be a NS LMXRB and 

* Based on observations with INTEGRAL, an ESA project with in- is one of the six cuiTently known Galactic Z so urces^ It dis- 

struments and science data centre funded by ESA member states (es- plays a complete Z pattern in the CC and HIDs (iKuuIkers et alJ 

pecially the PI countries: Denmark, France, Germany, Italy, Spain, and 

Switzerland), Czech Republic and Poland, and with the participation ' The other known Galactic Z sources are GX340+0, GX349+2, 

of Russia and the USA. Cyg X-2, GX 17+2 and Sco X-1 . 
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ll994tLTonkeret al.ll2002l) and shows secular shifts of the "Z". 
As most of the Z sources, it is located near the Galactic Centre 
which has rendered difficult its study due to source confusion 
and optical obscuration. Throughout the paper we assumed 
a distance of 8k pc and column density Nh = 3x10^^ cm~^ 
jAsai et al. 1994). Like all Z sources, GX 5-1 is a radio source 
(e.g. fender & Hendry. .2000.) with the radio emission most 
likely originating in a compact jet. The determination of the 
radio counterpart allows for extremely accurate position mea- 
surements. These have led to a most likely candidate for an 
infra-red companion ( Jonker et al. 2000). In the X-ray domain, 
up to now nei ther pulsations nor bursts have been detected (e.g. 
IVaughanetai.>1994) . 

Previous GX 5-1 hard X-ray data were almost always con- 
taminated by the nearby (40') BHC LMXRB GRS 1758-258. 
The first mission to clearly resolve GX5-1 from GRS 1758-258 
in the high energy domain was GRANAT (Paul et al. 1991). 
Observations from ART-P (~5-20keV), on board GRANAT, 
showed that GX 5-1 was ~30-50 times brighter than GRS 
1758-258 below 20keV. At higher energies SIGMA de- 
tected only GRS 1758-258 iSunvaev et al.ll99lHGilfanov et alJ 
Il993l) . 

In this paper we report the results of one year of monitoring 
of GX 5-1 with the INT Ernational Gamma-R ay Astrophysics 
Laboratory, INTEGRAL dWinkler et alJl2003l) . INTEGRAL is 
an ideal mission to study the hard X-ray emission of sources 
in the crowded Ga lactic Centre region. INTEGRAL/lBlS 
dUbertini et alJ2003 h has high sensitivity, about ^10 times bet- 
ter than SIGMA, coupled to imaging capability with 12' an- 
gular resolution above 20keV . In addition, two X-ray mon- 
itors, JEM-Xl and JEM-X2 jLund et alJl2003l) . provide the 
soft-X ray simultaneous coverage (3-35 keV, 3' angular reso- 
lution). Effectiveness of INTEGRAL in the study of Galactic 
bulge bright LMXRBs has been demonstrated in IPaizis et alJ 
C2003..2004.) using only initial data. 



2. Observations and data analysis 

INTEGRAL performed a Galactic Plane Scan (GPS) about ev- 
ery 12 days and Galactic Centre Deep Exposures (GCDEs) 
according to the Galactic Centre visibility. Among these ob- 
servation programs, GX 5-1 has been continuously in the 
INTEGRAL field of view around two distinct periods, April 
2003 (INTEGRAL Julian Date^ IJD~ 1200, hereafter data set 

1, -90ksec) and October 2003 (IJD~ 1380, hereafter data set 

2, ^77ksec). 

For the analysis we used the available X-ray monitor, JEM- 
X2 (hereafter JEM-X), for soft photo ns (3-35 keV) and the low 
energy IBIS detector, ISGRI (Leb run et al.l 120031) for harder 
photons (20-200 keV). We have not used PICsIT, the hard pho- 
ton IBIS detector dPi Cocco et al.l2003l) . as its peak sensitivity 
is above 200 keV, well beyond the 2 keV detection limit given 
by SIGMA. The spectrometer SPI dVedrenne etalJ '2003) has 
not been used as the presence of the nearby BH GRS 1758-258 



^ The INTEGRAL Julian Date is defined as the fractional number of 
days since January 1, 2000 (TT). IJD=MJD-51544. 



could contaminate the data. The separation of the two sources, 
40', is below the angular resolution of the instrument, ~ 2.5°. 

JEM-X and IBIS/ISGRI data were reduced using the Off- 
line Scientific Analysis (OSA) package, OSA 4.0. The descrip- 
tion of the algorithms used in the scientific analysis can be 
found in Westergaard et al., (2003) and Goldwurm et a l. ( 20031) 
for JEM-X and IBIS/ISGRI respectively. XSPEC 11.3.0 was 
used for the spectral analysis. The Comptonisation model by 
Nishimura et al. (1986) (COMPBB) was updated so that the 
cutoff energy of the model, originally at 70keV, would suit 
ISGRI spectra (i.e. moved up to 300 keV). 

We used JEM-X to build the hardness intensity diagram 
(HID) of GX 5-1. This was done using all the JEM-X data 
available with off-axis angle < 5° which led to 113 sci- 
ence windows (hereafter ScWs) around 2000 sec exposure 
each^. We extracted individual images from each ScW in 3 
en ergy bands (3-5 , 5-12 and 12-20 keV). A mosaicking tool 
bv lChenevez et al. (2004) was used to obtain the final JEM-X 
mosaic image. Within each ScW we extracted 100 s bin light- 
curves and used one ScW and 100 s data bins to build the HID. 
For the spectral fitting of JEM-X data it is preferable to have 
the source within 3° off-axis (better vignetting correction and 
signal-to-noise ratio). This led to 44 JEM-X ScWs. Spectra 
were grouped so that each new energy bin contained at least 
200 counts. Based on Crab calibration results, 3% systemat- 
ics have been applied between 5 and 22 ke V were the spectral 
analysis has been performed. 

In the case of ISGRI we selected the data within the totally 
coded field of view (FOV, <4.5°, to avoid the partially to to- 
tally coded FOV edge), resulting in a total of 90 ScWs (44 of 
which are in common with JEM-X). We extracted individual 
ScW images in 10 energy ranges with the boundaries 20, 30, 
40, 50, 60, 80, 100, 150, 300, 500 and lOOOkeV. The fluxes 
extracted in each image and energy band were used to build 
the light-curves of GX 5-1. Images from each ScW were com- 
bined in a final mosaic from which an overall ISGRI spectrum 
was extracted. Spectra were also extracted from each ScW with 
the Least Square Method and were further grouped so that each 
new energy bin had a minimum of 20 counts. A 5% systemat- 
ics has been applied to all channels. For the spectral fitting we 
used the OSA 4.0 ancillary response file (ARE) and redistribu- 
tion matrix (RMF) re -binned to 16 spectral channels between 
20 and 120keV. 

Simultaneous ISGRI and JEM-X fitting was carried out for the 
44 JEM-X ScWs. The remaining 46 ISGRI ScWs were not 
analysed separately because there was not enough statistics to 
study ISGRI spectra alone on a ScW basis. In the JEM-X - 
ISGRI simultaneous fit a cross-calibration factor was let free to 
vary. 



A science window is the basic unit of an INTEGRAL observation. 
Under normal operations one ScW will correspond to one pointing or 
one slew. To build the HID, we considered 113 pointings, 13 from the 
GPS (~2200 sec each) and 100 from the GCDE (-1800 sec each). 
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Fig. 1. Mosaic images combining the April 2003 and October 2003 data sets. Left panel: JEM-X mosaic image of GX 5-1 in 
the 5-lOkeV band (44 ScWs, ~76ksec). Right panel: IBIS/ISGRI mosaic image of G X5-1 in the 20-30 keV band (90 ScWs, 
~167ksec). The source IGR J17597-2201 has been recently discovered by INTEGRAL (Lutovinov et alJ2003l) . 
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Fig. 2. IBIS/ISGRI light curves of GX 5-1 in the 20-30 keV for the April 2003 and October 2003 coverage (left and rigth panel, 
respectively). Each point is one ScW (about 2000 s). Error bars are at la. The observed countrate corresponds to a variability 
between about 30 and 130mCrab. In the softer band (i.e. in JEM-X) the source is much brighter, around 1 Crab. 



3. Results 

3. 1. Overview of the data 

Figure [1] shows the JEM-X, 5-lOkeV, and IBIS/ISGRI, 20- 
30keV, mosaic image of the region of GX 5-1 (April 2003 
and October 2003 data sets combined). The nearby BHC 
GRS 1758-258 is rather dim in the softer energy range (JEM- 
X) and becomes much brighter in the harder range (ISGRI) 
where it is clearly disentangled from GX 5-1 . In the ISGRI mo- 
saic, we detect GX 5-1 with about 4 counts/sec and a detection 
significance of about 100. In the same image, GRS 1758-258 
is detected with an average countrate of about 2 counts/sec and 
a detection significance of about 50. For comparison, the av- 
erage count rate obtained from Crab, in different positions of 
the totally coded field of view, is about 75 counts/sec in the 
20-30 keV. In Fig. H the overall ISGRI 20-30 keV fight curve 
of GX 5-1 is shown. Each data point corresponds to one ScW 
in which GX 5-1 has been automatically detected by the soft- 



ware. This happens in 66 ScWs (out of 90) in the 20-30 keV 
band (with detected countrate between 2 and 10 counts/sec, i.e. 
between about 30 and 130mCrab) and in 10 ScWs in the 30- 
40keV band (less than 2 counts/sec, '^60mCrab). In the re- 
maining ScWs, GX 5-1 was not reaching the detection signif- 
icance of 3, set as a threshold for the automatic detection in a 
single ScW. The combination of all the ScWs (^2 ksec each) in 
a single mosaic (167 ksec, Fig.Q allows to reveal a hard X-ray 
emission above ~20keV. 

In the softer energy range covered by JEM-X, GX 5-1 is much 
brighter (around 1 Crab) and is detected in each 100 sec bin 
with an average of 27 counts/sec (3-5 keV), 56 counts/sec (5- 
12keV) and 5 counts/sec (12-20 keV). 



4 



Paizis et al.: Resolving the hard X-ray emission of GX 5-1 with INTEGRAL 




1.0 r 1 

40 50 80 1D0 120 140 

COUNTS/SEC 3 - 12 KEV 

Fig. 3. HID of GX5-1 with JEM-X (data set 1 ). The soft colour 
is the (5-12keV) to (3-5 keV) flux ratio. The graph was built 
using 100 s time bins that were smoothed with a three bin box- 
car: this means that each 100 s data value is replaced by the 
mean of the three neighbouring values. The smoothing was 
chosen for visualisation purposes as it makes the "Z" more ev- 
ident. The horizontal branch (upper part, soft colour > 2) and 
normal branch (lower part, soft colour < 2) are clearly visible. 

3.2. The hardness intensity diagrams 

In order to build the HID of GX 5-1 we used the JEM-X soft 
colour, defined as the (5-12 keV) to (3-5 keV) flux ratio, plot- 
ted versus the JEM-X intensity (defined as the 3-12keV flux). 

Figure |3] shows the HID we obtained for the first data set 
(April 2003). The graph was built starting from the 100 s bins 
that were smoothed to a 5 minute final bin for visualisation 
purposes. The horizontal branch, HB, and normal branch, NB, 
are quite evident while the flaring branch, FB, the last part of 
the "Z", is not visible. The FB could be hidden in the data 
since in GX 5-1 it is short and/or spread out by our choice 
of energy boundaries. Note the similarity with the HID de- 
rived from GINGAI ASM data jyan der Klis et alJ [l99l and 
Mir/Kvant TTM data (' Blom etalJl993l) . 

Figure |3 shows the "raw", i.e. non-smoothed, HIDs ob- 
tained for the two data sets (April 2003 and October 2003 
respectively). In the first data set the HB-NB vertex was set 
to (118,2.15) judging from the overall distribution of the data 
points. Then the slopes of the "Z" were calculated so that an 
equal number of data points were left on either side of the line. 
The deduced slopes are shown in Fig. |3 with a solid line. As 
the first data set has a cleaner "Z" pattern, we used this one 
as the "true" one to which the second "Z" was forced to bend. 
This was done shifting the HB-NB vertex of the second data 
set until the slopes of HB and NB became equal to those of the 
first data set. This led to a vertex in the second set at (115.1, 
2.4) which resulted in a clear secular shift from 2.15 to 2.4 in 
the soft colour of the vertex. Since the current data have an al- 
most continuous coverage from the HB through the NB while 
coverage of the FB (if any) is very poor, in this paper we focus 
on file HB and NB. 

In order to investigate the coiTelation between the source 
position in the "Z" pattern of the HID and the spectral be- 



haviour, we introduced a one dimensional parameter that mea- 
sures the position on the "Z": (^Dieters & van der Klis 2000[ 
and references therein). Strong evidence shows that the mass 
accretion rate M of individual Z sources increases from the 
top-le ft to the bottom-right of the "Z" pattern (|Hasinger e t alJ 
Il990l) . i.e. along the HB, NB and FB. Hence, the S z parameter 
is supposed to be monotonically related to the mass accretion 
rate. 

For each JEM-X data point (100 s bin and one ScW bin) in 
the HIDs, the two projections onto both branches and the dis- 
tance to each branch were computed. The branch at the shortest 
distance was chosen and after confirmation via visual inspec- 
tion, to ensure that data points did not jump from branch to 
branch outside the vertex area, the value of Sz was assigned. 
The scale was set so that the HB-NB vertex was at Sz = 1 
while Sz = and S^ = 2 were assigned at a 3-12 keV in- 
tensity of 40 counts/s (beginning of HB and end of NB respec- 
tively). In this frame, Sz<l means the source is on the HB 
(lower M) and l<Sz<2 means that the source is on the NB 
(higher A/). 

3.3. Energy Spectral Modeling 

We averaged all the ISGRI and JEM-X data to achieve the 
best statistics and tried to find the most physically reasonable 
spectral model to fit the entire data set. The ISGRI (average) 
spectrum was obtained from the mosaic image shown in Fig.[l] 
right panel, and the JEM-X total spectrum by averaging the 44 
JEM-X ScW spectra. In this case, the JEM-X spectrum was 
further grouped to have a number of energy bins comparable to 
ISGRI. 

The X-ray spectra of bright LMXRBs hosting a NS are gen- 
erally described as the sum of a soft and a hard component. Two 
different models are often adopted to describe t his composite 
emission. They are the so- called eastern mo del jMitsuda et alJ 
1984) and western model jWhite et alJl986l) . 

In the eastern model the softer part of the spectrum is a 
multi-colour disc describing the emission from the optically- 
thick, geometrically-thin accretion disc (XSPEC DISKBB 
model. Mitsuda et'al]ll984l) . The temperature at the inner disc 
radius kTi„ and the normalisation, linked to the inner disc ra- 
dius itself, are the parameters of the model. The harder part 
of the spectrum is a higher temperature Comptonised black- 
body describing the emission fro m the NS boundary layer 
(COMPBB model, Nishimura et an il986l) . Photons from the 
NS surface at a temperature kT;,b are up-scattered by a plasma 
of temperature kTg and optical depth r. In this case the nor- 
malisation is linked to the seed photon emitting area. 

In the western model a single temperature blackbody (from 
the optically thick boundary layer) describes the soft part of the 
spectrum (BB model in XSPEC). Comptonisation of soft seed 
photons (kTo) in the innermost region of the accretion disc by a 
corona of given temperature kTg and optical depth r describes 
the hard part (Titarchuk 1994, e.g. COMPTT model). 

We used the eastern and western model to fit the average 
JEM-X and ISGRI spectra simultaneously, and compared the 
results. In both cases a galactic absorption by a fixed hydrogen 
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Fig, 4, HID of GX5-1 with JEM-X for the two data sets. The scattered data points are 100 s bins while the points with error bars 
are one ScW bins. The average statistical error of the 100 s bins is shown in the bottom right of the graphs. The ScW bin values 
are computed as average of 100 s bins and the associated error is not the statistical one but a measure of how much the source 
moved along the "Z" track within the ScW. The solid line shows the deduced "Z" for the two data sets. Left panel: data set 1 
(April 2003). Right panel: data set 2 (October 2003). The "Z" shape is less evident in this data set with respect to the left panel. 
This change could be intrinsic to the source as the exposure time in the two data sets is comparable (58 ScWs in data set 1 and 
55 ScWs in data set 2). 
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Fig. 5. Best fit of the average ISGRI and JEM-X spectra us- 
ing the eastern model. The parameters of the fit are given in 
Tabled Residuals in terms of a and the single spectral compo- 
nents are shown. 

equiv alent column density of Nh = 3x 10^^ cm^^ dAsai et all 
11994 was added. The results of the fit are given in TableQland 
Fig. Eland li 

Both models can describe the spectral shape below 
~30keV equally well. However, in the western model, there 
is a significant residual in the ISGRI data above ^30keV, that 
is reasonably well described in the eastern model. This is be- 
cause the western model predicts an exponential cut-off above 
~10keV, while the eastern model describes the spectral "flat- 
tening" above ~20keV as a Comptonised hard-tail emission. 
We are not able to constrain the Comptonising plasma temper- 
ature (and consequently the cut-off energy) from our average 
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Fig. 6. Best fit of the average ISGRI and JEM-X spectra using 
the western model. Residuals in terms of a and the single spec- 
tral components are shown. A hard excess above ^30keV is 
visible. The parameters of the fit are given in Table [l] 



spectrum. In the eastern model, a good fit can be obtained with 
a plasma temperature up to about 50 keV. Either the cut-off is 
above the fitted energy range or our data are not of high enough 
quality to constrain the spectral break. We choose to fix the 
temperature to lOkeV as in this frame we can explain in a co- 
herent way the spectral evolution of GX 5-1 (see below). 

Adding a power-law component with fixed photon index 
r=2.5 in the western model results in a good fit of the ob- 
served hard-excess as shown in Fig.El In the 5-100 keV band 
the power-law component alone contributes to about 1.5% of 
the total luminosity. 

Figure|8]shows the photon spectrum obtained with the east- 
ern model parameters shown in Table [l] Since the eastern 
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Table 1. Best fit parameters for the average ISGRI and JEM-X 
spectra of GX 5-1. F5_20kcV and F2o-iookeV are the unab- 
sorbed fluxes in the 5-20 keV and 20-100 keV range respec- 
tively, in units of ergs^^cm^^; d.o.f. = degrees of freedom. 
The indicated errors are at la. No error means the parameter 
was fixed to the indicated value. The cross-calibration factor 
was frozen to 0.9 for ISGRl with respect to JEM-X. 





Eastern model 


Western model 


kTin or kT (keV) 


1.4 


1.7 


DISKBB or BB norm 


288 


0.1 


kTbb or kTo 


1.93±0.01keV 


0.87±0.02keV 


kXe 


lOkeV 


3.02±0.02 keV 


r 


0.37±0.01 


4.15±0.06keV 


COMPBB/COMPTT norm 


136 


1.6 


Red. 


1.12 (14 d.o.f.) 


6.21(13 d.o.f.) 


F5-20keV 


1.54x10"** 


1.54x10"** 


F20-100keV 


3.41x10"^" 


3.10x10"^" 




-5 10 20 50 100 

Energy [KeV] 

Fig. 7. Photon spectrum of GX 5-1 obtained with the west- 
ern plus additional power-law (r=2.5) model (Red. x^=1.2, 
12 d.o.f.). Dash: NS boundary layer emission (BB). Dot- 
dash: Comptonisation of innermost region of the accretion disc 
(COMPTT). Dot: additional power-law. Solid: total spectrum. 

model provides a physical interpretation for the hard tail, we 
focus on this model for the analysis and discussion of the spec- 
tral changes of GX 5-1. 

3.4. Spectral variability along the "Z" 

We studied the spectral variation of GX 5-1 along the HID us- 
ing the Sz parameter within the frame of the eastern model 
(with one ISGRl - JEM-X combined spectrum per ScW for a 
total of 44 ScWs). The soft component (DISKBB) was frozen 
to the values obtained in the best fit given in Table 1, for all 
the ScWs. Similarly, the plasma temperature kT^ was fixed at 
a value of lOkeV, since this parameter is strongly correlated 
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Fig. 8. Photon spectrum of GX 5-1 obtained with the eastern 
model described in Tabled Dash: soft, disc blackbody compo- 
nent. Dot-dash: hard, NS surface blackbody component with- 
out Comptonisation (r=0). A hard excess is evident. Dot (ba- 
sically overlapping with the solid line): NS blackbody compo- 
nent to which Comptonisation from a lOkeV, t=0.37 plasma 
was added. Solid: total spectrum (disc component and NS 
Comptonised component). The latter two are basically overlap- 
ping above lOkeV i.e. almost all the emission is coming from 
Comptonisation of the NS surface photons. 

with r. The normalisation of the hard component, kT;,;, and r 
were let free to vary from one Sew fit to the next as well as 
the instrument cross-calibration factor. An average value of re- 
duced x^ — 1.2 was obtained. We found a clear relation between 
the Sz value and the deduced kTbb for all the ScWs. To see if 
we could describe the spectral evolution of GX 5-1 in a more 
constraining way, we proceeded fitting the data letting in one 
case kTfib vary (with r fixed to 0.4) and in the other r vary 
(kThf, fixed to 2.0 keV). In both cases the normalisation was let 
free as well. Fixing kTf,f, resulted in a much worse fit (average 
reduced x^— 2.5) while fixing r gave an equally acceptable fit 
(average reduced x^—^-^)- This is consistent with the fact that 
kThb has a clear relation with the position of GX 5-1 along the 
"Z" and should not be frozen. Conversely, t shows a random 
distribution of values around 0.4 (extreme cases being 0. 1 and 
0.8), it is not the main parameter driving the spectral evolution 
along the "Z" and hence can be fixed. 

The combined JEM-X and ISGRI spectral fitting for two 
ScWs, one for the HB and the other one for the NB, are shown 
in Fig. 121 left and right respectively. The spectrum extracted 
from the HB is significantly harder than the NB one. This is 
true also for the other spectra and the hard X-ray emission is 
indeed characteristic of the HB. In fact we obtain that the hard 
X-ray emission (ISGRI countrate) is systematically linked to 
the position of the source in the CC diagram: Fig. 1101 shows 
the correlation among ISGRI 20-40 keV flux and JEM-X soft 
colour for the two datasets (off-axis angle <5°). The different 
symbols used refer to the position of the source in the "Z" track: 
the cases that have been identified with the HB have a higher 
ISGRI flux whereas the NB cases have a lower ISGRI flux. The 
second dataset, Fig.[lO|right panel, has more uncertain cases as 
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Fig. 9. ISGRI and JEM-X spectra for GX 5-1 in two different ScWs. The best fit and residuals in terms of cr are shown. The two 
components of the model (DISKBB and COMPBB) are also shown (dashed lines). The Comptonising component is frozen to 
kTe=10keV and r=0.4. Left panel: ScW 012000950010 (IJD=1377.57, UTC=2003-10-09). For this ScW Sz=0.50, i.e. GX 5-1 
is in the HB. kTb;,=2.35keV, red. x^=1.4 (112 d.o.f.). (The residuals obtained in the two ISGRI bins between 20-30 keV are not 
significant as they are due to known ISGRI systematics). Right panel: ScW 012200110010 (IJD=1381.63, UTC=2003-10-13). 
S2=1.56 i.e. GX 5-1 is in the NB. The spectrum is softer, only one ISGRI (> 20keV) data point met the requirement of a 
minimum of 20 counts per bin. kTbb=1.49keV, red. x^=1.5 (89 d.o.f.) 
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Fig. 10. Correlation among ISGRI 20-40 keV flux and JEM-X soft colour, intensity ratio (5-12keV)/(3-5keV), for the two 
datasets (April 2003 on the left and October 2003 on the right). Each point is one ScW. The different symbols used refer to the 
position of the source in the "Z" track: □ = clearly HB, = uncertain HB, A = HB or NB, x = clearly NB, * = uncertain NB, -|- 
= NB or FB. The big boxes indicate transitions between HB and NB (vertex S = 1). The two outliers in the left panel (assigned 
HB with low ISGRI intensity) are possibly due to systematics. 



expected given that its "Z" shape is less clean than for the first 
dataset. 

Figure ^3 shows the variation of the emitting hot NS sur- 
face area and temperature, kTbf,, along the "Z" (left and right 
panel respectively). While GX 5-1 moves from the HB to the 
NB (left to right in Fig. II 1> . the size of the emitting surface 
increases while its temperature decreases. 

4. Discussion 

Our observations constitute the first opportunity in which the 
Z source GX 5-1 could be monitored in a long term period in 



the hard X-ray domain without contamination from the nearby 
BHC GRS 1758-258. 

The source was thoroughly covered by INTEGRAL mainly 
during two periods (April and October 2003). The HIDs 
of GX 5-1 for the two periods show a clear secular shift. 
This is consistent with previous obser vations of GX 5-1 it- 
self as well as of other Z sources llKuulkers etalJll994t 
iKuuIkers & van der KlisI 1 19951 and references therein) for 
which a secular shift was detected. In our data, GX 5-1 cov- 
ers widely the HB and NB while very poor coverage of the 
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Fig. 11. Relations among the seed photons' properties and position of the source in the "Z". The electron plasma temperature kTe 
and the optical depth were frozen to lOkeV and 0.4 respectively. In both panels the two ScWs of Fig.|9lare highlighted. Left panel: 
relation between the blackbody emitting surface and Sz- The surface is in units of km^ and is computed as Trxd^ x normalisation, 
where d is the distance of the source in units of lOkpc, assumed to be 0.8. Right panel: relation between the derived blackbody 
temperature kT^f, and Sz- 



FB (if any) is seen, similarly to previous observations (e.g., 
iKuulkers etalJl994 . 

4. 1. Hard X-ray Emission 

The ISGRI and JEM-X average spectra show that GX 5-1 has 
a clear hard X-ray emission above 20keV. Assuming a dis- 
tance of GX 5-1 of 8kpc, we obtain luminosities of the order 
of Li_20koV--10^^ergs-i and L20-200koV'-10^^ ergs"!, 
fully consistent with v alues obtained for other Z sources 
llDi Salvo & Stellall2002h . With such luminosities, similarly to 
the case of the Z sources GX 11+2 and GX 349H-2, GX 5-1 
lies outside the so-called burster box (Barret etal. 2000), in 
a region that was originally believed to be populated by BH 
LMXRBs. Like BHs, NS LMXRBs can have a hard X-ray 
emission when they are bright. 

At least in the cases in which a cut-off is observed, the hard 
X-ray emission can be interpreted as thermal Comptonisation 
of soft photons in a hot, optically thin plasma in the vicinity 
of the neutron star. In the case of GX 5-1 we are not able to 
constrain the Comptonising plasma temperature. This suggests 
that either our data are not of good enough quality to measure 
the cut-off or that the cut-off is beyond the fitted energy range. 
In any case, we do obtain a good fit with relatively low plasma 
temperatures (lOkeV) hence non-thermal processes, even if not 
ruled out, cannot be claimed. 

In cases where the cut-off has not been observed up to about 
100 keV or higher (Di Salvo et al, 2000, 2001), extremely high 
electron temperatures are required, which is unlikely, given that 
the bulk of the emission of Z sources is very soft. Non atten- 
uated power-law tails dominating the spectra at high energies 
may be produced by mildly relativistic {v/c ^ 0.1) outflows 
with flatter power laws corresponding to higher optical depth 

Obtained extrapolating our INTEGRAL best fit model down to 
IkeV. 



of the scattering medium and/or higher bulk electrons veloci- 
ties, i n a way that is similar to thermal Comptonisation jPsaltid 
I2OOII) . Alternatively, X-ray components could originate from 
the Comptonisation of se ed photons by the non th ermal high 
energy-electrons of a jet llDi Salvo & Stellall2002l and refer- 
ences therein) or via synchrotron self-Compton emission di- 
rectly by the jet ( Markoff & Nowak 2004) . 

4.2. Spectral Variation 

With INTEGRAL, we are able to study the spectral state of 
GX 5-1 within one ScW seeing the source variability above 
20keV that has never been studied before. On the other hand, 
the softer part of the spectrum, dominated by the disc emission, 
cannot be constrained and thus we fixed it to the value found 
from the JEM-X average spectrum. 

The spectral changes of GX 5-1 along the "Z" can be ex- 
plained by the smooth variation of the physical properties of 
the seed photons originating from the hot NS surface (heated 
by the boundary layer). This sce nario is different from the 
more complicated ones that e.g. iHasinger et alJ ( Il990l) and 
lHoshi&Mitsudal ( ll99ll) proposed to explain GINGA data of 
Cyg X-2 and GX 5-1 respectivelv. .Hasinger et al. (1990) stud- 
ied the spectral variations of Cyg X-2 using both the eastern 
and western model. The temperature and the luminosity of both 
components (soft and hard) were let free to vary in both mod- 
els and a rather complicate interplay of the parameters of both 
components was describing the spectral variations of the source 
along the "Z" pattern. Hoshi & M ijsud_a ( 1991) used the east- 
ern model for GX 5-1 and the spectral changes were also in 
this case due to the variation of different parameters: in the NB 
the inner disc radius temperature kTi„ and the NS boundary 
layer temperature kT^b remained constant and spectral shape 
changes were explained by an increase of the emission area 
(from lower left to upper right NB). While the spectrum be- 
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came harder, in the HB, an additional parameter, kTi„, was let 
free to reproduce the trend of the data. 

Unlike GINGA, with INTEGRAL data we cannot constrain 
the variability of the disc component and thus we had to fix 
it to the values found from the average spectra. The normali- 
sation of this component (modelled with DISKBB in XSPEC 
terminology) led to an inner disc radius of Rm(cos6')^/^ 
~14(cos6')^/^ km for an assumed source distance of 8kpc 
where 6 is the angle between the normal to the disc and 
the line of sight. In the inner region of the disc, a hot outer 
layer affects the emergent spectra by Comptonisation. This 
leads to an observed temperature Tcoi which is higher than 
the effective temperature Teff- As a consequence the ob- 
served inner disc radius is underestimated by a factor of 
(TrniJ Tpff)^- With this hard ening factor correction assumed to 
be 1.5 ('F.bisa waetal.ll99lL we obtained an inner disc of about 
Re//(cos6i)i/2~30(cos6i)i/2 km. 

Not being sensitive to changes in the softer part of the spec- 
trum is a limitation in our study as we do expect the soft part of 
the spectrum (the disc) to vary but, for the first time, we were 
able to study the variability of GX 5-1 in the less explored hard 
X-rays, above 20 keV. We see that the hard component is indeed 
chatacteristic of the HB, namely the spectrum is harder at lower 
accretion rates. This seem s to be the general trend in LMXRBs 
llBarret & Vedrenn^l 19941) . Atoll sources are thought to be at 
a lower level of accreting rate than Z sources and indeed they 
are generally harder. This trend continues also within a class of 
source s: the hard components observed in the Z so urces GX 
17H-2 llDi Salvo et al. 2000..). Cvg X -2 llDi Salvo et al.. .2002) 
and GX349-H2 dOi Salvo et al J200lh decrease in intensity from 
HB to NB and FB, i.e. for increasing mass accreti on rate. An 
exception to this trend was found in Sco X-1 ( D'Amico et alJ 
I2OOII) where the presence of the hard component was not re- 
lated to the position of the source in the "Z" track. Theoretical 
interpr etations for this general trend have already been dis- 
cussed jlnogamov & Sunvaevll999HPopham & Sunvaevl200l[ 
iKluzniak & WagoneJl98.'C 

IPopham & SunvaevI 1I2OO II) have computed the boundary 
layer structure and its evolution with mass accretion rate M. 
They found a significant dependence of the boundary layer size 
(both height and radial extension) with M. Matter falls on a por- 
tion of the NS surface that is touched by the boundary layer. 
This part of the NS surface is thermalised and emits blackbody 
emission. As M increases, the area covered by the boundary 
layer will increase and so will the NS surface responsible of 
the blackbody emission. This is the same trend we found in 
the case of GX 5-1 (Fig. ^2 left): S^, i.e. M, increases and 
the emitting NS surface increases. This trend is more clearly 
visible in the HB than in the NB. A possible explanation for 
this could be that along the NB (with increasing mass accre- 
tion rate), the accretion disc, i.e. the soft part of the spectrum 
that we are not able to constrain with INTEGRAL, is starting 
to change. Besides, in the NB M is not expected to vary much 
(as it is already close to MEdd) so this may be also (part of) the 
reason why there is no clear trend in the NB. 

The 5-100 keV observed luminosity reached a maximum 
variation of a factor of ~3, unabsorbed Ls-iookcV— (0.6- 



1.7)xl0^*ergs"l (^). This is small compared to the NS 
emitting surface variation (factor of ~20, see Fig. ^2 left). 
Ls-iookeV is the total luminosity obtained from the compo- 
sition of soft and hard component between 5-100 keV but the 
variation in Ls-mokeV is to be attributed to the changes in the 
hot NS surface properties alone as the boundary layer domi- 
nates the emission and the disc as well as the parameters of 
the Comptonising plasma were kept fixed. In our interpreta- 
tion, any variation in Ls-iookoV is a consequence of changes 
in the properties of the hot NS surface that we modelled with a 
blackbody emission. Its luminosity is proportional to the emit- 
ting surface and temperature, Lbb ocR^xTf,f,^. So, if the lu- 
minosity stays nearly constant (factor of ~3) while the emit- 
ting surface increases (factor of '^20), then the temperature of 
the emitting surface will decrease. Thus, the anti-correlation 
that we found between Sz and the blackbody temperature kT^f, 
(Fig. ^2 right) may be qualitatively understood. 

5. Conclusions 

We have studied the X-ray emission of the Z source GX 5-1 
with INTEGRAL. A clear hard emission above ~20keV is de- 
tected and it may have the same origin of the hard components 
observed in other Z sources as it shares the property of being 
characteristic to the HB. We used the so-called "eastern" and 
"western" model to fit the average spectra of GX 5-1 . The east- 
ern model describes the spectral "flattening" above ^20keV as 
a Comptonised hard-tail emission whereas the western one pre- 
dicts an exponential cut-off above ~10keV and an additional 
power-law is needed to account for the hard X-ray emission. 
Since the eastern model provides a physical interpretation for 
the hard tail, we focussed on this model for the analysis of the 
spectral changes of GX 5-1. 

We interpret the spectral changes of GX 5-1 along the 
"Z" pattern in the hardness intensity diagram in terms of 
Comptonisation of varying soft photons (^2keV) by a hot 
plasma (lOkeV). The soft photons are interpreted as black- 
body emission from the part of the NS surface that is heated 
by the boundary layer The Comptonising plasma is the bound- 
ary layer optically thin plasma. When GX 5-1 moves down- 
wards in the "Z", the temperature and optical depth of the 
Comptonising plasma are not seen to change in the data. WTiat 
changes along the "Z" is the temperature of the optically thick 
emission from the hot NS surface that shows a steady decrease 
with increasing mass accretion rate. This may be a consequence 
of the gradual expansion of the boundary layer that we detect 
in the data. This trend is in agreement with theoretical studies 
( Po pham & Sunvaevll200ll) that indeed predict the expansion 
of the boundary layer surface with increasing mass accretion 
rate. 

With the INTEGRAL long term monitoring we will also 
have higher chances to catch GX 5-1 in the FB. The be- 
haviour of Z sources on the FB, believed to correspond to 
super-Eddington accretion, is quite complex and it will be in- 
teresting to see if it can be smoothly included in the HB and 

^ This range includes the 10% uncertainty due to the current cali- 
bration status. 
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NB scenario described here. Furthermore we will search for 
hard tails extending above 100 keV, in GX 5-1 as well as in the 
other Z sources and NS binaries of our monitoring program. 

The synergy of a long term monitoring and a multi- 
wavelength study should be the best way to understand the 
physics of hard-X ray emission of LMXRBs. Our collaboration 
is moving in this direction: regular coordinated INTEGRAL- 
RXTE observations of a sample of AtoU and Z sources (includ- 
ing GX 5-1) are being performed and simultaneous radio and 
optical observations are attempted. 
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